Dynamics of Multi-Player Games 
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We analyze the dynamics of competitions with a large number of players. In our model, n players 
compete against each other and the winner is decided based on the standings: in each competition, 
the Tnth ranked player wins. We solve for the long time limit of the distribution of the number of 
wins for all n and m and find three different scenarios. When the best player wins, the standings are 
most competitive as there is one-tier with a clear differentiation between strong and weak players. 
When an intermediate player wins, the standings are two-tier with equally-strong players in the 
top tier and clearly-separated players in the lower tier. When the worst player wins, the standings 
are least competitive as there is one tier in which all of the players are equal. This behavior is 
understood via scaling analysis of the nonlinear evolution equations. 

PACS numbers: 87.23.Ge, 02.50.Ey, 05.40.-a, 89.65.Ef 



I. INTRODUCTION 

Interacting particle or agent-based techniques are a 
central method in the physics of complex systems. This 
methodology heavily relies on the dynamics of the agents 
or the interactions between the agents, as defined on a 
microscopic level \^. In this respect, this approach is 
orthogonal to the traditional game theoretic framework 
that is based on the global utility or function of the sys- 
tem, as defined on a macroscopic level Q]- 

Such physics-inspired approaches, where agents are 
treated as particles in a physical system, have recently 
led to quantitative predictions in a wide variety of so- 
cial and economic systems |j, |j, |^. Current areas of 
interest include the distribution of income and wealth 
[a, Hi 13 J opinion dynamics llOl IllL ll^ , the propaga- 
tion of innovation and ideas 1131 , and the emergence of 
social hierarchies 14, 15, 16lll7|. 

In the latter example, most relevant to this study, com- 
petition is the mechanism responsible for the emergence 
of disparate social classes in human and animal commu- 
nities. A recently introduced competition process p4Lll7l| 
is based on two-player competitions where the stronger 
player wins with a fixed probability and the weaker player 
wins with a smaller probability [ig. This theory has 
proved to be useful for understanding major team sports 
and for analysis of game results data ^ . 

In this study, we consider multi-player games and ad- 
dress the situation where the outcome of a game is com- 
pletely deterministic. In our model, a large number of 
players n participate in the game, and in each competi- 
tion, the TTith ranked player always wins. The number 
of wins measures the strength of a player. Furthermore, 
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the distribution of the number of wins characterizes the 
nature of the standings. We address the time-evolution 
of this distribution using the rate equation approach, and 
then, solve for the long-time asymptotic behavior using 
scaling techniques. 

Our main result is that there are three types of stand- 
ings. When the best player wins, m = I, there is a 
clear notion of player strength; the higher the ranking 
the larger the winning rate. When an intermediate player 
wins, 1 < m < n, the standings have two tiers. Players in 
the lower tier are well separated, but players in the upper- 
tier are all equally strong. When the weakest player wins, 
TO = n, the lower tier disappears and all of the players 
are equal in strength. In this sense, when the best player 
wins, the environment is most competitive, and when the 
worst player wins it is the least competitive. 

The rest of this paper is organized as follows. We in- 
troduce the model in section II. In Section III, we analyze 
in detail three-player competitions, addressing situations 
where the best, intermediate, and worst player wins, in 
order. We then consider games with an arbitrary number 
of players and pay special attention to the large-n limit 
in Section IV. We conclude in section V. 



II. THE MULTI-PLAYER MODEL 

Our system consists of N players that compete against 
each other. In each competition n players are randomly 
chosen from the total pool of players. The winner is 
decided based upon the ranking: the mth ranked player 
always wins the game [Fig. ^. Let ki be the number 
of wins of the ith ranked player in the competition, i.e., 
ki > ■ ■ ■ > km > ■ ■ • > fcn, then 

[ki, . . . , K„i, ■ ■ ■ Kn) ^ (rCi, . . . , Km -|- 1, . . . , Knj. (1 j 

Tie- breakers are decided by a coin-toss, i.e., when two or 
more players are tied, their relative ranking is determined 



FIG. 1: Illustration of the multiplayer game with n = 5 and 
m — 3. 



in a completely random fashion. Initially, players start 
with no wins, fc = 0. 

These competition rules are relevant in a wide variety 
of contexts. In sports competitions, the strongest player 
often emerges as the winner. In social contexts and es- 
pecially in politics, being a centrist often pays off, and 
furthermore, there are auctions where the second highest 
bidder wins. Finally, identifying wins with financial as- 
sets, the situation where the weakest player wins mimics 
a strong welfare system where the rich support the poor. 

We set the competition rate such that the number of 
competitions in a unit time equals the total number of 
players. Thence, each player participates in n games per 
unit time, and furthermore, the average number of wins 
(fc) simply equals time 



(fc) - t. 



(2) 



At large times, it is natural to analyze the winning rate, 
that is, the number of wins normalized by time, x — k/t. 
Similarly, from our definition of the competition rate, the 
average winning rate equals one 



(x) = 1. 



(3) 



Our goal is to characterize how the number of wins, or 
alternatively, the winning rate are distributed in the long 
time limit. We note that since the players are randomly 
chosen in each competition, the number of games played 
by a given player is a fluctuating quantity. Nevertheless, 
since this process is completely random, fluctuations in 
the number of games played by a given player scale as the 
square-root of time, and thus, these fluctuations become 
irrelevant in the long time limit. Also, we consider the 
thermodynamic limit, N ^ oo. 



III. THREE PLAYER GAMES 

We first analyze the three player case, n = 3, because it 
nicely demonstrates the full spectrum of possibilities. We 
detail the three scenarios where the best, intermediate, 
and worst, players win in order. 



A. Best player Avins 

Let us first analyze the case where the best player wins. 
That is, if the number of wins of the three players are 
fci > fc2 > fca, then the game outcome is as follows 



Let fk{t) be the probability distribution of players with 
fc > wins at time t. This distribution is properly nor- 
malized, J2k fk = 1; ^^^d it evolves according to the non- 
linear difference-differential equation 



(/fc-i-F'fc-i - fkFf,) 



(5) 



+ 



{fl^iFk-i 



flFk) 



+ {:){fL,-fk) 



dt [l 
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Here, we used the cumulative distributions Fk = J2j=o fj 
and Gk — X) i=fe+i fj of players with fitness smaller than 
and larger than fc, respectively. The two cumulative dis- 
tributions are of course related, Fk + Gk-i = 1- The first 
pair of terms accounts for games where it is unambigu- 
ous who the top player is. The next pair accounts for 
two-way ties for first, and the last pair for three way ties. 
Each pair of terms contains a gain term and a loss term 
that differ by a simple index shift. The binomial coeffi- 
cients account for the number of distinct ways there are to 
choose the players. For example, there are (j^) = 3 ways 
to choose the top player in the first case. This master 
equation should be solved subject to the initial condition 
/fc(0) = 5k fi and the boundary condition f-i{t) — 0. One 
can verify by summing the equations that the total prob- 
ability is conserved ^ J2k fk = 0, and that the average 
fitness (fc) = ^j, fc/fc evolves as in ^, d{k)/dt = 1. 

For theoretical analysis, it is convenient to study the 
cumulative distribution Fk ■ Summing the rate equations 
(O , we obtain closed equations for the cumulative distri- 
bution 



dFk 
dt 



^3{Fk - Fk-i)Fl^, 



(6) 



-3(Ffc - Fk-ifFk-i 



<Fk-Fk^if. 



Here, we used fk = -Ffc+i — Fk- This master equation is 
subject to the initial condition -F'fe(O) = 1 and the bound- 
ary condition F_i(t) = 0. 

We are interested in the long time limit. Since the 
number of wins is expected to grow linearly with time, 
fc ^ i, we may treat the number of wins as a continuous 
variable, Fk._^ ^ Fk - ^ + ^fp^ + ■■ ■ . Asymptoti- 



'k-i 
cally, since ^ ex i" 



and 



afc2 



oc t 



etc., second- and 



higher-order terms become negligible compared with the 
first order terms. To leading order, the cumulative dis- 
tribution obeys the following partial differential equation 



^ + 3F'^ 

dt + ^ dk 



0. 



(7) 



From dimensional analysis of this equation, we anticipate 
that the cumulative distribution obeys the scaling form 



(fcl,fc2,fc3) -> (fcl + I,fc2,fc3)- 



(4) 



Ffe(i)~$(fc/i) 



(8) 




FIG. 2: The stationary distribution of winning rates p2^ for 
the case n — 3, m = 1. 



with the boundary conditions $(0) — and $(00) — 1. 
In other words, instead of concentrating on the number 
of wins k, we focus on the winning rate x — k/t. In the 
long time Umit, the cumulative distribution of winning 
rates $(x) becomes stationary. Of course, the actual dis- 
tribution of winning rates (j){x) also becomes stationary, 
and it is related to the distribution of the number of wins 
by the scaling transformation 



A(t)~i"V(fcA) 



(9) 



with (t>{x) — $'(a;). Since the average winning rate equals 
one (Q , the distribution of winning rates must satisfy 



1 = 



dxx^'{x). 



(10) 



Substituting the definition into the master equation 
(0, the stationary distribution satisfies 



$'(x)[3$2-x] =0. 



(11) 



There are two solutions: (i) The constant solu- 
tion, $(a;) = const, and (ii) The algebraic solution 
^six) = {x/Sy^^. Invoking the boundary condition 
limaj^oo ^{x) = 1 we find [Fig. E] 



^{x) 



(a;/3)i/2 a;<3 
1 a; > 3. 



(12) 



One can verify that this stationary distribution satis- 
fies the constraint (|10|) so that the average winning rate 
equals one. This result generalizes the linear stationary 
distribution found for two player games .17] . 

Initially, all the players are identical, but by the ran- 
dom competition process, some players end up at the top 
of the standings and some at the bottom. This directly 
follows from the fact that the distribution of winning 



rates is nontrivial. Also, since ^(a;) 



"^/^ as a: ^ 0, the 



distribution of winning-rate is nonuniform and there are 
many more players with very low winning rates. When 
the number of players is finite, a clear ranking emerges. 




t 



FIG. 3: Schematic illustration of the number of wins k ver- 
sus time t when the best player wins. Shown are results of 
simulations with 20 players. 



and every player wins at a different rate. Moreover, after 
a transient regime, the rankings do not change with time 
[Fig.El. 

We note that in our scaling analysis, situations where 
there is a two- or three-way tie for first do not contribute. 
This is the case because the number of wins grows linearly 
with time and therefore, the probability of finding two 
players with the same number of wins can be neglected. 
Such terms do affect how the distribution of the number 
of wins approaches a stationary form, but they do not 
affect the final form of the stationary distribution. 



B. Intermediate player wins 

Next, we address the case where the intermediate 
player wins, 

(fcl,fc2,fc3)->(fcl,fe + l,fc3). (13) 

Now, there are four terms in the master equation 
dfk /3\ /2^ 



dt 



ifk-lFk-lG, 



fc-1 



1/ VI 

(/, 

[fk-lPk-l - fk Pk 



fkFkGk) (14) 



A {fk-iGk-i - fi 



iGk) 



+ (3) ifk-i - fk) 



The first pair of terms accounts for situations where there 
are no ties and then the combinatorial prefactor is a 
product of the number of ways to choose the interme- 
diate player times the number of ways to choose the best 
player. The next two pairs of terms account for situa- 
tions where there is a two-way tie for best and worst, 
respectively. Again, the last pair of terms accounts for 
three-way ties. These equations conserve the total prob- 
ability, J2k fk = ^, and they are also consistent with ^. 



o 



1/4- 




FIG. 4: The stationary distribution of winning rates 1171 for 
n = 3, m — 2. 



Summing the rate equations (|14|l . we obtain closed 
equations for the cumulative distribution 



dFk 
dt 



-&{Fk - Fk-i)Fk-iGk-i 

-3{Fk — Fk-i) {Fk^i + Gk-i) 
-{Fk-Fk-if. 



(15) 



For clarity, we use both of the cumulative distributions, 
but note that this equation is definitely closed in Fk be- 
cause of the relation Gfc = 1 — Fk+i- Taking the con- 
tinuum limit and keeping only first-order derivatives, the 
cumulative distribution obeys the following partial differ- 
ential equation ^ -I- 6F(1 — -f )§f- = with the bound- 
ary conditions Fq — and linife^oo F^ = 1. Substituting 
the definition of the stationary distribution of winning 
rates (|S)| into this partial differential equation, we arrive 
at 



$'(a;)[6$(l-$)-2;] = 0, 



(16) 



an equation that is subject to the boundary conditions 
$(0) = and lima^^oo ^(x) — 1. There are two so- 
lutions: (i) The constant solution, ^{x) = const, and 
(ii) The root of the second-order polynomial ^s{x) = 

^(1 



2 y-.- -^/l — 2x/3). Invoking the boundary conditions, 
we conclude [Fig. ^ 




X < xo 

X > Xq. 



(17) 



As the nontrivial solution is bounded $s(x) < 1/2, the 
cumulative distribution must have a discontinuity. We 
have implicitly assumed that this discontinuity is located 
at xo < 3/2. 

The location of this discontinuity is dictated by the 
average number of wins constraint. Substituting the sta- 
tionary distribution H17|l into H1Q(I then 



1 = 



dxx^\x) +xo[l - ^(xo)]. 




t 



FIG. 5: Schematic illustration of the number of wins k versus 
time t when the intermediate player wins. Shown are results 
of simulations with 20 players. 



In writing this equality, we utilized the fact that the sta- 
tionary distribution has a discontinuity at Xo and that 
the size of this discontinuity is 1 — <I>o- Integrating by 
parts, we obtain an implicit equation for the location of 
the discontinuity 



1 



Xo 



dx^{x). 



(18) 



Substituting the stationary solution H17() into this equa- 
tion and performing the integration, we find after several 
manipulations that the location of the singularity satis- 
fies the cubic equation Xq(^xq — |) =0. The location of 
the discontinuity is therefore 



9 

Xo = 7T- 



(19) 



This completes the solution (|17|) for the scaling function. 
The size of the discontinuity follows from $0 = ^(a^o) = 
1/4. 

There is an alternative way to find the location of the 
discontinuity. Let us transform the integration over x 
into an integration over $ using the equality 



dx$(a;)-t- / d$a;($). (20) 



a;o$o 



This transforms the equation for the location of the dis- 
continuity (|18|l into an equation for the size of the jump 



a;o(l-$o) 



d^x{^). 



(21) 



<&) we arrive at the cubic 

1 = 6$o-9$^+4«>i^- The 

from which we conclude 



Substituting x{^) = 6$(1 - 
equation for the variable $0, 
relevant solution is $0 — 1, 
Xo = 9/8. For three-player games, there is no partic- 
ular advantage for either of the two approaches: both 
(|18|l and H21|) involve cubic polynomials. However, in 
general, the latter approach is superior because it does 
not require an explicit solution for ^(x). 



The scaling function corresponding to the win-number 
distribution is therefore 



<"=> = U'-r 



-1/2 



-da; 

4 1 8 



where d{x) denotes the Kronecker deha function. The 
win-number distribution contains two components. The 
first is a nontrivial distribution of players with winning 
rate x < xq and the second reflects that a finite fraction 
of the players have the maximal winning rate x = xq. 
Thus, the standings have a two-tier structure. Players in 
the lower tier have different strengths and there is a clear 
differentiation among them [Fig.|H]. Players in the upper- 
tier are essentially equal in strength as they all win with 
the same rate. A fraction <i>o = -j belongs to the lower 
tier and a complementary fraction 1 — $0 = | belongs to 
the upper tier. Interestingly, the upper-tier has the form 
of a condensate. We note that a condensate, located at 
the bottom, rather than at the top as is the case here, 
was found in the diversity model in Ref. jl7'|. 



C. Worst player ■wins 

Last, we address the case where the worst player wins 

MM 



{ki,k2,k'i) -^ {ki,k2,k3 + 1). 



(22) 



Here, the distribution of the number of wins evolves ac- 
cording to 



dfk 
dt 



{fk-iGk-i - fkGf.) 
(/fe-iGfe-i - fk Gk) 
) (./ti - fl) ■ 



(23) 



This equation is obtained from |(S| simply by replacing 
the cumulative distribution Fk with Gk ■ The closed equa- 
tion for the cumulative distribution is now 



dF\ 
dt 



-3{Fk - Fk-i)Gk_i 
— 3(Ffe — Fk-i) Gk~i 



(24) 



-{Fk-Fk-if. 

In the continuum limit, this equation becomes 
^ + 3(1 — F)"^^ = 0, and consequently, the stationary 
distribution satisfies 

$'(x)[3(l-$)^-a;] =0. (25) 



Now, there is only one solution, the constant ^{x) = 
const, and because of the boundary conditions <I>(0) = 
and \\ra.x^oo^(x) = 1, the stationary distribution is a 
step function: <I>(a;) = 1 for x > Xq and ^{x) = for 



<|) 



FIG. 6: The stationary distribution of winning rates 1261 for 
the case n = m = i. 



X < Xq. In other words, $(a;) = Q{x — xq). Substitut- 
ing this form into the condition (|1U|I . the location of the 
discontinuity is simply xq — 1, and therefore [Fig. |^ 



$(x) = e(x-l) 



(26) 



where Q{x) is the Heaviside step function. When the 
worst player wins, the standings no longer contain a 
lower-tier: they consist only of an upper-tier where all 
players have the same winning rate, (l){x) = S{x — 1). 



IV. ARBITRARY NUMBER OF PLAYERS 

Let us now consider the most general case where there 
are n players and the ?7ith ranked player wins as in (Q . It 
is straightforward to generalize the rate equations for the 
cumulative distribution. Repeating the scaling analysis 
above, Eqs. Hll|l and (|16|l for the stationary distribution 
© generalize as follows: 



$'(a;)[C$"-™(l-$)' 



0. 



(27) 



The constant C equals the number of ways to choose the 
TOth ranked player times the number of ways to choose 
the 771 — 1 higher ranked players 



G = 



n — 1 

TTl — 1 



{n — m)\{m — 1)! 



(28) 



Again, there are two solutions: (i) The constant solu- 
tion, $'(x) — 0, and (ii) The root of the {n — l)th-order 
polynomial 



C^' 



'{l-^y 



(29) 



We now analyze the three cases where the best, an inter- 
mediate, and the worst player win, in order. 
Best player wins (m = 1): In this case, the stationary 
distribution can be calculated analytically. 



$(x) 



(a;/n)i/("-i) x < n] 



X > n. 



(30) 



One can verify that this solution is consistent with Q. 
We see that in general, when the best player wins there is 
no discontinuity and $o = 1- As for three-player games, 
the standings consist of a single tier where some players 
rank high and some rank low. Also, the winning rate 
of the top players equals the number of players, a;o — 
n. In general, the distribution of the number of wins is 
algebraic. 

Intermediate player wins (1 < m < n): Based on the 
behavior for three player games, we expect 



$(a;) 



1 X > Xq. 



(31) 



Here, ^s{x) is the solution of H29|l . Numerical simula- 
tions confirm this behavior [Fig. |2|. Thus, we conclude 
that in general, there are two tiers. In the upper tier, all 
players have the same winning rate, while in the lower 
tier different players win at different rates. Generally, a 
finite fraction $o belongs to the lower tier and the com- 
plementary fraction 1 — <I>o belongs to the upper tier. 

Our Monte Carlo simulations are performed by simply 
mimicking the competition process. The system consists 
of a large number of players N, all starting with no wins. 
In each elemental step, n players are chosen and ranked 
and the mth ranked player is awarded a win (tied players 
are ranked in a random fashion). Time is augmented by 
1/A^ after each such step. This elemental step is then 
repeated. 

The parameters xq and $0 characterize two important 
properties: the maximal winning rate and the size of each 
tier. Thus, we focus on the behavior of these two parame- 
ters and pay special attention to the largc-n limit. Substi- 
tuting the stationary distribution H31() into the constraint 
(|10|l . the maximal winning rate xq follows from the very 
same Eq. (|18|) . Similarly, the size of the lower tier follows 
from Eq. (|21|l . In this case, the latter is a polynomial of 
degree n + I, so numerically, one solves first for $0 and 
then uses (|29|l to obtain xq. We verified these theoretical 
predictions for the cases n = 4 and n = 10 using Monte 
Carlo simulations [Fig. [7] . 

For completeness, we mention that it is possible to 
rewrite Eq. (|21|) in a compact form. Using the definition 
of the Beta function 



Jo 



d$$'^ 



$ 



,m — 1 



B{n-m + l,m) (32) 
(n — m)\{m — 1)! 



= c- 



we relate the definite integral above with the combinato- 
rial constant in H28|l . Substituting the governing equation 
for the stationary distribution (|29|l into the equation for 
the size of the lower-tier (|21|l gives 



c-i = $^'-" 



(l-*o)^ 



d$$" 



'(1-$)" 



(33) 




FIG. 7: The stationary distribution of winning rates "l>(x) for 
n — 4 (top) and n = 10 (bottom). Shown are Monte Carlo 
simulation results with A'' = lO'' particles at time t = 10^. 
The circles are the theoretical predictions for the maximal 
winning rate xo and the size of the lower tier "l>o. 



Using the relation (|32|l . we arrive at a convenient equa- 
tion for the size of the lower tier $0 



fi$$"-™(l_$) 



ni—1 ^^"" 



'{i-<^oy 



(34) 



This is a polynomial of degree n + 1. 

Let us consider the limit n —^ 00 and to — > 00 with 
the ratio a = m/n kept constant. For example, the case 
a = 1/2 corresponds to the situation where the median 
player is the winner. To solve the governing equation for 
the stationary distribution in the large-n limit, we esti- 
mate the combinatorial constant C using Eq. H28|l and 
the Stirling formula n! ~ {2Tmy/^{n/ey\ Eq. ^ be- 
comes 



na 



27r(l-a) \l-a 



$ 



1-$ 



~ X. 



(35) 



Taking the power 1/n on both sides of this equation, 
and then the limit n ^ 00, we arrive at the very simple 
equation. 



$ 



1 



1 -$ 



= 1. 



By inspection, the solution is constant, <I> = 1 
$0 = 1 ~ CK and employing the condition (x) 



(36) 

a. Using 
= 1 yields 



1-a 



^ 



1/a 



FIG. 8: The infinite-n limit. From Eq. I|87fl . the points 
{xo, $o) all lie on the curve $ = {x — l)/x. 



the location of the condensate 



Xo = 1/a, 



$n 



1 



(37) 



This result is consistent with the expected behaviors 
Xq — > oo as a — > and xo(a = 1) = 1 (see the worst 
player wins discussion below). Therefore, the station- 
ary distribution contains two steps when the number of 
players participating in each game diverges [Fig. |S] 



<^{x) 



X <0 
a < X < a~^ 



(38) 



a 



< X. 



The stationary distribution corresponding to the num- 
ber of wins therefore consists of two delta-functions: 
(f>ix) = (1 — a)S{x) + a6{x — 1/a). Thus, as the number 
of players participating in a game grows, the winning 
rate of players in the lower tier diminishes, and eventu- 
ally, they become indistinguishable. 

For example, for n — 10, the quantity $o is roughly 
linear in a and the maximal winning rate xq is roughly 
proportional to a~^ [Fig. [7]. Nevertheless, for moder- 
ate n there are still significant deviations from the limit- 
ing asymptotic behavior. A refined asymptotic analysis 
shows that $0 — (1 ~ <^) ~ 1/0(1 — a) In n/n and that 
Xo — (1 — $o)~^ |2y|- Therefore, the convergence is slow 
and nonuniform (i.e., a-dependent). Despite the slow 
convergence, the infinite-n limit is very instructive as it 
shows that the structure of the lower-tier becomes trivial 
as the number of players in a game becomes very large. 
It also shows that the size of the jump becomes propor- 
tional to the rank of the winning player. 

It is also possible to analytically obtain the station- 
ary distribution in the limit of small winning rates, 
X — > 0. Since the cumulative distribution is small, $ — > 0, 
the governing equation (|29|l can be approximated by 
(j^n-m — j,_ As a result, the cumulative distribution 
vanishes algebraically 



$(x) -- X"-™, 

as X — > 0. This behavior holds as long as m < n. 



(39) 



Worst player wins {m — n): In this case, the roots of 
the polynomial (|29|l are not physical because they cor- 
respond to either monotonically increasing solutions or 
they are larger than unity. Thus, the only solution is 
a constant and following the same reasoning as above 
we conclude that the stationary distribution is the step 
function (|26|l . Again, the upper tier disappears and all 
players have the same winning rate. In other words, there 
is very strong parity. 

We note that while the winning rate of all players ap- 
proaches the same value, there are still small differences 
between players. Based on the behavior for two-player 
games, we expect that the distribution of the number of 
wins follows a traveling wave form Fk{t) -^ U{k — t) as 
t — > cx) Q. As the differences among the players are 
small, the ranking continually evolves with time. Such 
analysis is beyond the scope of the approach above. Nev- 
ertheless, the dependence on the number of players may 
be quite interesting. 

Let us imagine that wins represent wealth. Then, the 
strong players are the rich and the the weak players are 
the poor. Competitions in which the weakest player wins 
mimic a strong welfare mechanism where the poor ben- 
efits from interactions with the rich. In such a scenario, 
social inequalities are small. 



V. CONCLUSIONS 

In conclusion, we have studied multi-player games 
where the winner is decided deterministically based upon 
the ranking. We focused on the long time limit where sit- 
uations with two or more tied players are generally irrele- 
vant. We analyzed the stationary distribution of winning 
rates using scaling analysis of the nonlinear master equa- 
tions. 

The shape of the stationary distribution refiects three 
qualitatively different types of behavior. When the best 
player wins, there are clear differences between the play- 
ers as they advance at different rates. When an interme- 
diate player wins, the standings are organized into two 
tiers. The upper tier has the form of a condensate with 
all of the top players winning at the same rate; in con- 
trast, the lower tier players win at different rates. In- 
terestingly, the same qualitative behavior emerges when 
the second player wins as when the second to last player 
wins. When the worst player wins, all of the players are 
equal in strength. 

The behavior in the limit of an infinite number of play- 
ers greatly simplifies. In this limit, the change from upper 
tier only standings to lower tier only standings occurs in a 
continuous fashion. Moreover, the size of the upper tier is 
simply proportional to the rank of the winner while the 
maximal winning rate is inversely proportional to this 
parameter. 

In the context of sports competitions, these results are 
consistent with our intuition. We view standings that 
clearly differentiate the players as a competitive environ- 
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merit. Then, having the best player win results in the 
most competitive environment, while having the worst 
player win leads to the least competitive environment. 
As the rank of the winning player is varied from best to 
worst, the environment is gradually changed from highly 
competitive to non-competitive. This is the case because 
the size of the competitive tier decreases as the strength 
of the winning player declines. 

In the context of social dynamics, these results have 
very clear implications: they suggest that a welfare strat- 
egy that aims to eliminate social hierarchies must be 
based on supporting the very poor as all players become 
equal when the weakest benefits from competitions. 

Our asymptotic analysis focuses on the most basic 
characteristic, the winning rate. However, there are in- 
teresting questions that may be asked when tiers of equal- 



strength players emerge. For example, the structure of 
the upper tier can be further explored by characterizing 
relative fluctuations in the strengths of the top players. 
Similarly, the dynamical evolution of the ranking when 
all players are equally strong may be interesting as well. 
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